Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN M?eX style file vl.4) 



Statistical characteristics of the observed metal systems 
and problems of reionization 

M. Demianski^'^, A.G. Doroshkevich^, 

llnstitute of Theoretical Physics, University of Warsaw, 00-681 Warsaw, Poland 
2Department of Astronomy, Williams College, Williamstown, MA 01267, USA 

3Astro Space Center of Lebedev Physical Institute of Russian Academy of Sciences, 117997 Moscow, Russia 



Accepted Received in original form ... . 



ABSTRACT 

We analyze basic properties of about 200 metal systems observed in high resolution 
spectra of 9 quasars by Boksenberg, Sargent & Ranch (2003). The measured Doppler 
parameters for the hydrogen and carbon lines are found to be different by a factor 
~ ^ mc/iTiH what indicates the domination of the thermal broadening and high 
degree of relaxation for majority of the observed metal systems. The analysis of the 
mean separation of the metal systems confirms that they can be located in regions with 
the typical proper size Rd ^ 0.1 — 0.25h~^ Mpc what corresponds to the typical sizes 
of protogalaxies with the baryonic masses Mf, w 10^ — 10^'^ Mq. The metal abundances 
observed at redshifts z = 3 — 6 within both IGM and galaxies are found to be quite 
similar to each other what indicates the strong interaction of young galaxies with the 
IGM and the important role of dwarf satellites of host galaxies in such interaction. 
On the other hand, this metal abundance can be considered as the integral measure 
of nuclear reactions within any stars what in turn restricts the contribution of the 
stars to the creation of ionizing UV background. Available now estimates of these 
abundances demonstrate that attempts to explain the reionization of the Universe as 
a byproduct of the process of synthesis of metals only are problematic and a significant 
contribution of unobserved and/or non thermal UV sources seems to be required. 

Key words: cosmology: reionization of the Universe — quasars: absorption: metal 
abundance. 



1 INTRODUCTION 

One of the most promising way to study the earlier stages 
of the process of galaxy formation and reionization of the 
Universe is the analysis of the Ly-a forest and metal sys- 
tems observed in spectra of the farthest quasars. The great 
potential of such approach was discussed already by Oort 
(1981, 1984) just after Sargent et al. (1980) established the 
intergalactic nature of the forest. Indeed, the majority of HI 
absorbers are associated with the small scale distribution of 
intergalactic matter (IGM) along the line of sight at red- 
shifts z > 2, when the IGM is not yet strongly clustered and 
its observed characteristics can be more easily interpreted. 
In turn, the observed metal systems are evidently associ- 
ated with galaxies and they characterize the rate of galaxy 
formation and their environment. The available Keck and 
VLT high resolution observations of the absorption lines pro- 
vide a reasonable database and allow one to apply statistical 
methods for their analysis. This analysis can be also supple- 
mented by comparison of basic properties of metal systems 



with those available for galaxies at corresponding redshifts. 
These problems were discussed in many papers (see, e.g., 
Bergeron et al. 1992; Lanzetta et al. 1995; Tytler 1995; Le 
Brune et al. 1996) and more recently in Boksenberg, Sargent 
& Rauch (2003); Adelberger et al. (2005b) and Scannapieco 
et al. (2006). 

The enrichment of the IGM by metals can be related 
to the starburst driven outflows (Pettini et al. 2001; Frye, 
Broadhurst & Benitez 2002) from the main galaxy or its 
dwarf satellites and to the stripping of these satellites. Both 
these processes occur but their efficiency is different. Thus, 
rare richer multicomponent systems are usually associated 
with the periphery of massive galaxies (see, e.g., Adel- 
berger et al. 2005b; Songaila 2006). Their complex structure 
demonstrates that they were formed by a multi step process 
that involved both explosive and stripping metal injections. 
On the contrary, the commonly encountered poorer metal 
systems that are related to extended low density clouds are 
more probably linked with the dwarf satellites of the central 



© 0000 RAS 



2 Demianski & Doroshkevich 



galcixy. Detailed discussion of this problem can be found in 
Scannapieco et al. (2006). 

The comparison of Doppler parameters measured for 
the HI and metals in the same absorption systems allows 
also to compare the influence of the IGM temperature and 
the macroscopic motions on the structure of the clouds, to 
evaluate the degree of relaxation of compressed matter and 
to select the small fraction of clouds with the domination of 
supersonic macroscopic motions. 

The separate problem is the process of reionization of 
the Universe. Observations of the farthest quasars (Becker 
et al. 2001; Djorgovski et al. 2001; Fan et al. 2004) show that 
the reionization of the IGM had just boon completed at 2: ~ 
6 — 7. On the other hand, the measurements of temperature 
fluctuations of the CMB by WMAP (Hinshaw et al. 2006) 
suggest that the reionization of the IGM has began at least 
at 2 ~ 10 — 12. Simulations indicate that the first stars and 
supernovae can bo formed already at z ~ 40 — 50 (Reed ct 
al. 2005) but for the standard cosmological model at 2: > 10 
the fraction of matter accumulated by high density objects 
("galaxies") is negligible (< 10~®) and therefore reionization 
earlier than at z > 10 requires some exotic sources of the 
UV radiation such as, for example, antimatter or unstable 
particles (see, e.g, Cohen, De Rujula, Glashow 1998; Bambi 
& Dolgov 2007; Preese et al. 2007). 

It is commonly believed that the general features of the 
reionization process are already established with a reason- 
able reliability (see, e.g., Tumlinson et al. 2004; Madau 2007; 
Schaerer 2007). At z < 3.5 — 4 the observed radiation of 
quasars dominates the UV ionizing background, but it pro- 
vides not more than 20 30% of the UV background at 
z > 5. Now it is expected that the UV radiation required for 
the reionization can bo produced by the joint action of galax- 
ies dominated by Population III and Population II stars (see, 
e.g., Tumlinson et al. 2004; Scheierer 2007) and by sources of 
non thermal radiation such as AGNs, miniquasars or black 
holes (see, e.g., Madau & Rees 2001; Meiksin 2005). How- 
ever, the relative contribution of stars, gamma ray bursts 
(GRBs) and/or non thermal sources of UV radiation has not 
been determined yet. Now the main discussed problems are 
focused on indirect estimates and observational restrictions 
of efficiency of such sources (see, e.g., Dijkstra, Haiman & 
Loeb 2004; Meiksin 2005; Choudy et al. 2007) and the analy- 
sis of their interaction with the environment (see, e.g., Iliev, 
Shapiro & Raga 2004). 

However we still have not seen neither the Pop HI stars 
nor the first galaxies dominated by such stars and, so, the 
actual properties of such objects remain unknown. It seems 
that all galaxies observed at high redshifts -including both 
the Ly-breack galaxies (LBGs) and Ly-a emitters (LAEs) 
- are dominated by the Pop II stars, and only small part of 
LAE demonstrate some features expected for the first galax- 
ies. Because of this it can be expected that the domination 
of Pop III stars could take place during only a short pe- 
riod of evolution of the first galaxies and these stars acts 
as a trigger for the faster formation of Pop II stars and for 
the transformation of the first galaxies into the LAEs and 
LBGs (Schneider et al. 2002; Ricotti & Ostriker 2004; Smith 
& Sigurdsson 2007; Tornatore et al. 2007; Karlsson ct al. 
2007). Thus the actual properties of the first galaxies could 
be quite similar to the observed properties of the LBGs and 
LAEs and the contribution of stars (and particularly of the 



Pop III stars) to the production of the ionizing UV back- 
ground and the process of reionization can be quite limited 
(see, e.g., ShuU & Venkatesan 2007) while the contribution 
of alternative non thermal sources of UV luminosity can be 
essential. 

Special problem is the existence of significant popula- 
tion of LAEs at redshifts z < 5 — 6. It is evident that the life 
time of envelopes of HI surrounding these galaxies is finite 
owing to the combine action of the observed outer UV back- 
ground and internal sources of UV radiation. For the same 
reasons the recent formation or successive reconstruction of 
such envelops from the highly ionized IGM seem to be quite 
problematic. On the other hand, if LAEs represent special 
population formed from low ionized IGM at high redshifts 
then it is necessary to explain the unexpectedly slow evo- 
lution of survived objects. In this case a clear redshift vari- 
ations of the number density of LAEs and their observed 
properties should be observed. 

Here wc analyze the sample of 908 high resolution CIV 
absorbers presented in Boksenberg, Sargent & Ranch (2003) 
and compare the measured properties of CIV and HI lines. 
This approach allows us to demonstrate that the majority 
of metal systems are associated with relaxed clouds. Our in- 
direct estimates of the sizes of metal systems are consistent 
with direct measurements of Adelberger et al. (2005b) and 
Scannapieco et al. (2006). Our estimates of the CIV abun- 
dance in the IGM are similar to that observed within LBGs 
and LAEs, what demonstrates the important contribution 
of dwarf galaxies in the reionization. At the same time, the 
estimates of the metal abundance show that the radiation 
of stars only cannot reionize the Universe and the contribu- 
tions of these stars and non thermal sources can be at least 
comparable. 

In this paper we consider the spatially flat ACDM model 
of the Universe with the Hubble parameter and mean density 
given by: 

i/'(«) = H^Qm{l + zf[inA/flm{l + 2)"^1 , 

{nb{z)} = 2.4 • 10"'^(1 + zf{nbh^ /0.02)cm~^ , (1) 

{Pm{z)} = ^^rn{l + ^)^ Ho = 100ft km/s/Mpc . 

Here Qm ~ 0.3&iflA ~ 0.7 are the dimensionless density of 
matter and dark energy, ~ 0.02 and h — 0.7 are the di- 
mensionless mean density of baryons, and the Hubble con- 
stant. For z > 1 the influence of A-term in (1) becomes 
negligible and we will write the Hubble parameter as 

H{z) w HoVa^il + zf^ . (2) 
2 THE DATABASE 

The present analysis is based on 9 high resolution spec- 
tra listed in Boksenberg, Sargent & Rough (2003). The full 
sample of metal systems contains 908 CIV lines majority 
of which represent the internal structure of rich metal sys- 
tems. From these spectra we select 193 metal systems at 
2 < z < 4.4. For 93 of these systems we have also parame- 
ters of the HI absorbers. 
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Figure 1. For the sample of 93 metal systems with Nh > 
functions bcH = \J mc /mH bmet/t>H, and Nh are plotted vs. 



As is seen from Fig. 1 the hydrogen and carbon col- 
umn densities are weakly correlated and even for richer hy- 
drogen absorbers with Nhi > 10^® cm~^ there are metal 



systems with 2 • ^ 



< Nc < 6 ■ lO^cm"^ . 154 



systems with Nhi < 10 crn^ and Nciv < 
can be related to the intergalactic ones while 32 systems 
with Nhi > lO^^cm"^ are probably linked with intervening 
galaxies. The separation of these metal systems can be com- 
pared with results obtained for 268 publicly available metal 
systems observed with intermediate resolution. 

As is seen from Fig. 2, the redshift distribution of the 
metal systems is non homogeneous and the majority of sys- 
tems are concentrated at 2 < z < 3.5 . This means that some 
of the discussed here characteristics are derived mainly from 
this range of redshifts. At 2 > 3.5 the statistics of lines is 
not sufficient. Detailed discussion of the observed character- 
istics of these metal systems can be found in Boksenberg, 
Sargent & Rough (2003). 



3 OBSERVED CHARACTERISTICS OF 
METAL SYSTEMS 

The sample of 154 metal systems with Nc < 10^*cm~^, 
and Nh < lO^^cm"^ can be used to characterize the metal 
systems actually related with the IGM. For this sample the 
redshift variations of the four mean observed characteris- 
tics, namely, the column density of CIV, (Nc), the Doppler 
parameter, (bmet), the mean abundance of carbon in the 
IGM, (f2c), and the mean separation of metal systems, 
{D,ep} = {Dl^pZl^ ■'"'), Zi = (l)/4 are plotted in Fig. 2 
for 2 < 2: < 4 . Here the mean abundance of the carbon in 
IGM is estimated as follows. 
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Figure 2. For the sample of 154 metal systems the redshift distri- 
bution of the metal systems, fmet, the mean Doppler parameter. 



(bmet), the mean carbon abundance in the IGM, {CIq) = {flc^l) 
and the mean system separation, {D*^^) = (Dgepzl ) are plot- 
ted vs. redshift z. 



{^c) = 



rricNc 8tvG 



The redshift variations of these parameters are fitted by 
expressions 



(bmet) = (11.2±l)km/s, 

(Nc) = (Nczi) = 3 • 10'^(1 ± 0.3)cm-^ 



^4 = (l)/4, 



(n*c) = {^czl) = [4 ± 2]/i-^ • 10" V^-^/O-S , (3) 

(DLp) = (i).ep2i-^') = (27.4±5.8)/i-V0-3/^^- Mpc, 

respectively. These estimates are consistent with those ob- 
tained in Boksenberg et al. (2003) for 'complex systems'. 
They are also quite similar to the estimates of Scannapieco 
et al. (2006) and Songaila (2001). For all 193 systems the 

mean separation decreases to 

(D,%) = {Dsepzl^'') = (23.8 ± 3.4)/r 1 ^JoI/Q^Mpc . (4) 

For 32 systems with Nhi > lO^*" cm" which can be re- 
lated to intervening galaxies the hydrogen and metals abun- 
dance increase up to: 



(Ig Nc /Nhi) = -4.2 ± 0.6, (Ig A^c) = 14 ± 0.4 , 



(5) 



what indicates the low concentration of carbon as compared 
with its solar abundance {\gNc /Nhi)q ~ —3.5 (Allende 
Prieto et al. 2002). The mean separation of these objects is 



(Utep) = (DsepZi^'') = (93 ± 24)/i"V0-3/nmMpc , 

what substantially exceeds the previous ones. For 26 Ly- 
limit systems with Nhi > lQ^''cm~'^ this separation in- 
creases up to 

(L>:ep) = {Dsepzl '"') = (111 ± 31)/i-V0-3/nmMpc , (6) 
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Figure 3. For the sample of 193 metal systems PDFs for the 
Doppler parameter, P{xfi), the carbon column density, P{xc), 
and abundance, P{xn), and the separation of metal systems, 
P{xsep), are plotted together with fits (7). 

and the carbon abundance drops down to 

{Ig Nc/Nhi) = -4.6 ±0.2. 

For the 268 metal systems selected from spectra observed 
with intermediate resolution the mean separation is similar 
to (6). 

The observed probability distribution functions, PDFs, 
for the Doppler parameter, P{b), the carbon column density, 
P(Nczl), the carbon abundance, P{Q,c), and the separation 
of metal systems, P{Dsepz\'''^), are plotted in Fig. 3. Such 
choice of variables allows us to suppress the artificial rcdshift 
variations of the PDFs caused by the redshift evolution of 
the mean characteristics of these systems. These PDFs are 
well fitted by the following functions 

Pfrt(xc) ~ exp[-(xc - 0.8)Vl-2], xc = \gNc/{lgNc) , 

PfitiXb) W exp(-1.85a;6), Xb = bmet/{bmet) > 1 , 

-P/it(a;n) w 3.5exp(-2a;n), xn = zl / {^c zl) , (7) 
Pfitixa) w 1.9exp(-a;s), Xs = Dsepzl '^^ /{Dsepzl '^^) , 

where again Z4 = (l)/4. 

For the sample of 193 metal systems the main fraction of 
observed systems (~ 70%) is linked with rich absorbers with 
\g{NHi) = 15± 1 and only 10% of systems can be related to 
the Lyman-limit systems with \g{NHi) > 17. At the same 
time, 75% of absorbers with \g{NHi) > 14 do not contain 
metals above the measured limit lg{Nc) « 12. These results 
confirm that the observed pollution of the IGM is caused by 
rare sources of nictals. 

It is interesting also to compare the separations of metal 
systems with different richness. Thus for subsamples of 91 
poorer systems with Nc < lO^^cm"^ and 102 richer systems 



with Nc > 10^"^ cm'^ the mean separations are quite simi- 
lar, {Dsepzl'"^^) ~ 32ft^^Mpc. This result is well consistent 
with the mean separation obtained for the full sample of 193 
systems (4). It emphasises similarity in the spatial distribu- 
tions of poorer and richer systems concentrated within rare 
clouds. 

In spite of limited range of observed redshifts and the 
large scatter of estimates (3), growth with time of (fic) and 
(Nc) points out to the progressive enrichment of IGM by 
metals. These results arc consistent with those obtained by 
Boksenberg et al. (2003) and Scannapieco (2006) and can 
be related to the evolution of galaxies at the same redshifts. 
However, Songaila (2001) observed very slow redshift vari- 
ations of the carbon abundance flc — (5.2 ± 1.7) ■ 10^** at 
2 < z < 5.5 while Schaye et al. (2003) discussed the slow evo- 
lution of unexpectedly high carbon abundance flc « 2-10"'^. 

At the same time, the weak redshift variations of both 
(Nc) and Nc/Nhi (5) indicate that for the stronger systems 
related to the intervening galaxies the carbon abundance is 
redshift independent. The properties of metal systems are 
very important for reconstruction of galaxy evolution at high 
redshifts and they deserve more detailed observational inves- 
tigation. 



4 DISCUSSION 

Results presented in previous Section allow us to test the 
degree of relaxation of matter compressed within such ab- 
sorbers and to obtain indirect estimates of spatial character- 
istics of metal systems. Finally we can estimate the possible 
contribution of stars and the process of metal production to 
the radiation that reionized the Universe. 



4.1 Degree of absorbers relaxation 

Comparison of the Doppler parameters of metal systems, 
bmet, and accompanied hydrogen lines, 6h, allows to clarify 
the complex problem of the internal structure of absorbers. 
Thus, numerical simulations show that absorbers can be re- 
laxed along one or two shorter axes and the width of the 
absorption lines depends upon the thermal broadening, 6r, 
the differential Hubble flow, bHu, and turbulent macroscopic 
velocities, bmt, 

— 1/2 

where br oc depends upon the mass of ion, m;, while 

two other components are the same for all ions. 

The relative contribution of these factors varies from 
absorber to absorber (see, e.g., Theuns, Schaye & Haehnelt 
2000; Schaye 2001). For majority of absorbers the possible 
contribution of Hubble flow is naturally linked with the re- 
tained expansion of clouds and is more essential for weaker 
absorbers. Its contribution depends upon the (unknown) rel- 
ative orientation of absorber expansion and the line of sight. 
On the other hand, the supersonic macroscopic motions are 
rapidly transformed to the shock waves and so they disap- 
pear. Their life time is relatively short and existence over 
an extended period is quite problematic what implies that 
is a more probable result. 
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Comparison of the Doppler parameters of metal systems 
and accompanied hydrogen lines allows to discriminate be- 
tween contributions of thermal and various large scale bulk 
and turbulent motions. Indeed, macroscopic and turbulent 
velocities are the same for all ions and in that case we can 
expect that 



bn ~ be, bch = ^Jrnc /niHbc /bn ~ \/T2 ^ 3.5 . 



(8) 



On the contrary, for any given temperature of the gas the 
thermal velocities of ions depend upon their masses and in 
this case hch determined by (8) should be close to one. 

For 93 systems plotted in Fig 1 comparison of these 
Doppler parameters shows that 



{bch) = \/rac/mH{bc/bH) w 1-1 ± 0.5 , 



(9) 



and only for 8 absorbers 2 > bch > 1-5. This result indicates 
that the thermal velocities dominate within majority of ab- 
sorbers in the sample, br > bmt- As is seen from Fig. 1 a 
scatter of this ratio increases for small Nc < W^'^crn^^ what 
indicates smaller degree of relaxation for poorer metal sys- 
tems. This scatter increases if the pollution of IGM is caused 
by the ejection of metals from closest galaxies when the re- 
maining velocity of ejection increases the Doppler parame- 
ter of metals and makes it larger than that for the hydrogen 
component (see, e.g., Songaila 2006). 

These results imply the approximate hydrostatic equi- 
librium of compressed matter along the line of sight or at 
least along the shorter axis of clouds and agree well with 
inferences of CarswcU, Schayc & Kim (2002); Telfcr ct al. 
(2002); Simcoe, Sargent & Rauch (2002), (2004); Boksen- 
berg, Sargent and Rauch (2003); Bergeron & Herbert-Fort 
(2005). In particular, comparison of the Doppler parame- 
ters measured for HI, CIV and OVI (Carswell, Schaye & 
Kim 2002) verifies also that as a rule the macroscopic (tur- 
bulent) velocities are subsonic. 

These observational results strongly suggest that most 
absorbers are long-lived partly gravitationally bound and 
partly relaxed and they are composed of both DM and bary- 
onic components. 



4.2 Links between metal systems and galaxies 

Natural links between the metal systems and galaxies was 
confirmed by direct observations by Adelberger (2005a, b). 



where for majority of systems with Nc 



lO'^cnr^ the 



proper sizes of absorbing regions at 2; = 2 — 3 were estimated 
as Rabs ~ O.OSMpc for CIV Unes and as Rabs ~ 0.4Mpc for 
OVI lines. At the same time, Scannapieco et al. (2006) found 
that the metal systems are clustered within "bubbles" with 
the comoving size Rtt ^ 2/t^^Mpc. Analysis of the observed 
absorber separations allows us to test these results. 

The mean comoving separation of clouds is determined 

as 



{Dscp) — 



l + z 



3nciiz)Sciiz) 



(10) 



where Ud & Sd are the proper number density and sur- 
face area of the clouds projection on the sky. For the 
metal systems with Nhi > 10^® cm~^ we have (Dsep) « 
can be expected that such systems are 
related to gravitationally bounded intervening galaxies and 



for them Sd only weakly depends on redshift, Sd ~ const. 
Therefore, for such systems (ud) oc (1 -|- 2)^'^^ what implies 
that at redshifts 2; ~ 2 — 4 the comoving number density of 
galaxies increases with time as Ugai oc (l-t-a)""'^^. Compar- 
ison of the observed number density of galaxies in the SDSS 
and 2dF surveys (Doroshkevich et al. 2004) 

(11) 



^/i^Mpc" 



{ngai{z = 0)) ~ 10 
and of LBG galaxies at z 

{Ugaliz = 3)) Ah^ 



3 (Adelberger 2005a), 



10"^Mpc"^ , 



(12) 



indicates similar evolution, Ugai oc {l+z)~^^^ . Using (6& 12) 
we get for the proper mean surface area and size of such 
metal systems 



(Sci) ~ 0.05/1- 



■^Mpc^ 



(Rci) = \/5cz/7i-~ O.m-^Mpc. (13) 



This estimate is similar to the proper sizes of absorbing re- 
gions obtained by Adelberger (2005b) at 2 = 2 — 3 for CIV 

lines. 

For weaker metal systems we can expect progressive 
growth of the proper cloud surface area with time, (Sd) oc 
(1 -|- z)~^ and for /3 > 0.25 a successive decrease with time 
of their comoving number density (ud) oc (1 + z)^''^''^^ . 
This means that due to clouds expansion some fraction of 
the weak metal systems goes under the observational limit 
Nc < lO^^cm"^. Such effect was found for the Ly-a forest 
(Demiariski et al. 2006). 

Using the estimate (3) for the mean system separation 
we get for the mean proper surface area and size of such 
systems 

{Sd) « 0.2/1-^2- 



(J^c^> w 0.25/1-^2-"'" Mpc, (14) 

where 24 = (1 + 2)/4. This estimate of Rd is of about two 
- three times smaller than the values found by Adelberger 
et al. (2005b) for weak lines OVI and the comoving size of 
"bubbles", Rbb ~ 2h~^Mpc, found by Scannapieco et al. 
(2006) at redshifts 2 = 2-3. 

Formation of galaxies with the baryonic mass Mt im- 
plies accumulation of matter from the protogalaxy with the 
proper size 

1/3 



Rp^g « (6F/7r)^/^ « OAz^^Mpc 



Mb 0.02 
IOIOM0 Qbh? 



(15) 



what quite well agrees with (14 for Mb ~ 10^ M© and with 
observations of Adelberger et al. (2005b) and Scannapieco 
ct al. (2006) for Mb ~ 10^° M©. 

Now the pollution of IGM by metals is related either 
to ejection of metals from galaxies after explosions of super- 
novae or to the stripping of satellites of the central galaxy. It 
is most probable that both processes are equally important. 
Thus, as was noted in Sec. 3.2, the complex structure of 
richer systems implies the multistep metal ejection what is 
naturally connected with SN explosions in the large central 
galaxy. On the contrary, at large distances from the central 
galaxy ejection of metals from satellites and their stripping 
dominate. 

The pronounced correlation between the observed 
galaxies and metal systems reflects the expected interac- 
tion between the large and small scale perturbation when 
former ones modulate the process of galaxy formation and 
amplify the inhomogeneities in the spatial galaxy distribu- 
tion. At the same time, the noted above similarity in the 
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spatial distributions of richer and poorer metal systems in- 
dicates similarities in their genesis and their close links with 
the ordinar observed galaxies. 

4.3 Probable sources of reionization 

Measurements of the carbon abundance in the IGM at high 
rcdshifts allows us to clarify some aspects of the process of 
reionization of the Universe. Thus, the concentration of met- 
als in the older galaxies observed at redshifts z ^ 6 — 7 and 
higher is the integral measure of the contribution of nuclear 
reactions within any stars to the creation of the ionizing UV 
radiation. On the other hand, the comparison of the carbon 
abundances within gala:xies and in the IGM characterizes 
the efficiency of the processes of star explosion and ejection 
of metals from the host galaxies and their dwarf satellites. 
It is important that this approach estimates the integral ac- 
tion of all stellar sourscs ignoring the often discussed (see, 
e.g. Choudhury et al. 2007) but badly known details of the 
reionization process such as the mass function of the first 
galaxies and Pop. Ill stars, the rate of the star creation, rar 
diative feedback or the redshift variations of the transmitted 
flux and the electronic depth. 



Table 1. Galcixies at z > 3 



'z) 


M* 




Ref. 






lO^MoMpc-s 




- 3 


0.1 


0.3 


Lai et al. (2007) 


- 3 


1 


3. 


Lai et al. (2007) 


- 5 


1 


1.5 


Verma et al. (2007) 


- 6 


1 


1 - 7 


Yan et al. (2006) 


- 6 


10 


2.5 


Eylis et al. (2007) 


■5-6 


50 - 500 


8. 


Wiklind et al. (2007) 


- 7 


1-10 


1.6 


Labbe et al. (2006) 



M, and p* are the mean stellar mass and the stellar mass density 
for the observed sample of galaxies. 



(2007) (Table 1) which arc related to older massive com- 
pact galaxies with the comoving number density Ugai « 
1.4 • 10~®Mpc~^, the proper sizes Rg ~2kpc and ages 0.2- 
1 Gyr. It is expected that the majority of these stars were 
formed at z > 9 and they represent the objects actually 
responsible for the reionization. 

For this sample of galaxies we get for the stars and metal 
abundances 



4-3.1 Contribution of observed galaxies 

As is well known, the transformation of hydrogen to car- 
bon produces w 7.3MeV of energy per one baxyon what 
is equivalent to w 5 • 10^ of UV photons with energy 
Ery ss 13.6eV . The synthesis of O, N, Si and other metals 
produces the UV radiation with the similar efficiency. Ac- 
cording to more detailed estimates (see, e.g., Tumlinson et 
al. 2004) the value of ionizing photons per one baryon is 
smaller and varies with the mass of Pop. Ill stars in the 
range N-, ~ {2 - 8) ■ 10''. 

At the same time, the observed strong ionization of hy- 
drogen at 2 = 5 — 7 implies the generation already at such 
redshifts of at least one UV photon per baryon {Nt,^ > 1). 
Of course, this is the minimal value and in some papers (see, 
e.g. Dijkstra, Haiman, Loeb 2004; Madau 2007) production 
of extra (up to 10) UV photons per baryon is discussed. 
Thus, if the reionization was caused by the production of 
metals within stars then the minimal abundance of metals 
at 2 ~ 6 must be at least 



6.7-10" 



= Cl^Z « 1.3-10" 



/esc 



where iVt-y > 1 is the minimal number of UV photons per 
one baryon required for the reionization and fesc < 1 is 
the mean fraction of UV photons escaping from galaxies. 
The allowance for the complex spectral distribution of the 
generated UV photons, ionization of He and the heating and 
slow recombination of the IGM increases this estimate of the 
f2min by a factor of 2 3. 

Available now estimates of the stellar component in 
galaxy populations observed at redshifts z > 3 arc sum- 
marized in Table 1 . Here observations of Lai et al. (2007) 
relate to the two populations of LAEs (younger and older) 
and the scatter of measured p, is determined by variations 
of the selected samples of galaxies. 

The most interesting estimates for our goals came 
from recent observations of 11 galaxies by Wiklind et al. 



O.lZp 



(17) 



and fi, is significantly below the present day values 
Qi^rniz = 0) w 4 - 10"^ (Pukugita et al. 1998). It is evi- 
dent that this result is consistent with flmin (16) only for 
an unlikely value Nb-y/ fesc ~ 1- 

On the other hand, any such galaxy can ionize baryons 
in a volume with the mass Mi and comoving size Ri where 



Mi~20ee//M,, €eff 

P f 3Mi V^' / M, 0.02 



fesc Z 

■0X0.120 5-105 ' 

(. 



(18) 



1/3 



So the moderate size of ionized bubble corresponds to the 
small fraction of the ionized IGM, 



fi = MiUgal/iPb) ~ 10 €e// 



M, 0.02 

10" Mo Q.bh'^ ' 



(19) 



what illustrates the limited efficiency of ionization caused by 
the nuclear reactions within such galaxies. Of course, as was 
noted above, at redshift z = 3 the density of LBG galaxies 
(12) is larger by a factor of 10^ but for them M, < 10"Mq 
(Table 1) and finally we get again fi <C 1. 

Estimates of the star radiation (16, 19) include the fac- 
tor fesc- Its available estimates are quite uncertain as they 
strongly depend upon the internal structure of galaxy and 
the immediate environment of the UV sources. These esti- 
mates vary from fesc ~2 - 3% to fesc ~10% (Iliev et al. 
2004; Dijkstra et al. 2004; Meiksin 2005). For example, for 
the observed LAEs fesc and so majority of the UV radi- 
ation is absorbed within the galaxy. Thus, we can conclude 
that either all observations cited above and Eq. (17) strongly 
underestimate the metal production and the UV radiation 
of high redshifts galaxies or the reionization is related with 
some unobserved thermal and/or non thermal sources of UV 
radiation. 
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4-3.2 Possible contribution of dwarf galaxies 

The differences between estimates (16) and (17) will de- 
crease il the majority ol ionizing UV radiation is related 
with the dwarf unobserved galaxies. Such preferential format 
tion of dwarf galaxies is natural in the CDM cosmology and 
is consistent with observations (see, e.g., Yan & Windhorst 
2004; Yan ct al. 2006; Stark et al. 2007). It is also consis- 
tent with the shape of the UV luminosity function observed 
at redshifts « = 4 — 6 (Bouwens et al. 2007). Moreover, re- 
construction of the initial power spectrum in Demianski et 
al. (2006) indicates the possible excess of power at scales 
M < 10* M0 what can additionally stimulate earlier forma- 
tion of dwarf galaxies. 

The contribution of dwaxf galaxies is also restricted ow- 
ing to the observational limitations for the small fraction of 
dwarf galaxies survived at 2: = (sec also Adolborgor 2005a). 
Non the less, the important role of dwarf galaxies in the 
reionization is supported by the observed abundance of CIV 
in the IGM. All published estimates of carbon abundance 
at redshifts 3 < z < 5.5 are close to Qc ~ (5.2 ± 1.7) • 10~* 
(Songaila 2001; Pettini ot al. 2003; Boksenberg ct al. 2003; 
Scannapieco et al. 2006; Ryan- Weber et al. 2007) and are 
consistent with (3) in the range of scatter. The slow redshift 
evolution of the abundance implies that the main enrichment 
of the IGM takes place at high redshifts z > 6. Moreover, the 
carbon abundance in the IGM is also close to that observed 
in galaxies at the same redshifts (17) 

f2*C = ii*metZc /Zmet — 0.251, met ~ 2.6 ■ 10 ^ . (20) 

As was noted in the previous subsection these metals are 
presumably concentrated within bubbles with the comoving 
size ~ 2/i"^ Mpc. This size significantly exceeds the typi- 
cal size ~ 2 kpc observed for galaxies at tiic same redshifts 
(Wiklind et al. 2007) but it is close to the size of protogalaxy 
with the baryonic mass Mt ~ 10^° — 10^^ M©. In turn, such 
patchy spatial distribution of metals in the IGM implies that 
the process of metal enrichment is dominated by metal ejec- 
tion from the dwarf satellites of the host galaxy in the course 
of their hierarchical clustering. This inference is close to that 
formulated in Songaila (2006) where it was noted that the 
formation of weak CIV systems cannot be easily understood 
in terms of high velocity galactic wind. Observation of GRB 
at redshift z = 6.3 (Haislip et al. 2005) can be considered as 
a possible example of such processes. This discussion shows 
that the real metal abundance at > 6 exceeds the esti- 
mate (17) at least by a factor of 2. On the other hand, such 
scenario implies also the existence at z = of the invisible 
dwarf galaxies in the vicinity of massive galaxies of early 
types. 

4-3.3 Possible contribution of nonthermal sources 

This discussion shows that even taking into account all un- 
certainties of observations the attempts to explain the reion- 
ization as a consequence of nuclear reactions within stars 
seem to be highly problematic (see also Adelberger 2005a; 
Meiksin 2005). The non thermal sources of UV photons such 
as AGNs and black holes formed both in prcgalactic systems 
and within first galaxies produce ~ 50Mev per one accreted 
baryon what is 7 times more effective than the nuclear re- 
actions. This means that even accretion of a small fraction 



of baryons, 

nacr ^ 10-V/e.c ^ 1.4 • 10~''n,//e.c , (21) 

will produce the same effect as the observed stars (17). This 
verifies that such non thermal sources of UV radiation can 
be considered as very promising ones and that they can be 
actually responsible for the reionization (see, e.g., Madau 
&Rees 2001; Meiksin 2005; Reed 2005; Ciardi et al. 2006; 
Madau 2007). 

The observational restrictions of such non thermal 
sources of the UV radiation are poor. Thus, Meiksin notes 
that the comoving space density of galaxies with the black 
holes is ni,h ~ 2 • 10~'*Mpc~'' what is ^ 1% of galaxies at 
z = 0. However, this value is 10 times more than the comov- 
ing number density of early massive galaxies (Wiklind et al. 
2007) which demonstrate some features of the AGN activ- 
ity. A faint AGN was directly observed at z = 5.44 (Douglas 
ct al. 2007) . Restrictions of the contribution of such sources 
discussed in Dijkstra ct al. (2004) were criticized by Meiksin 
(2005). The possible contribution of exotic sources such as 
the unstable fraction of DM particles or antimatter can be 
noticeable at a > 20 - 30 (Cohen et al. 1998; Bambi & Dol- 
gov 2007; Freez at al. 2007) but it cannot dominate at the 
period of recombination. Its possible impact will be soon 
tested by the Planck mission. 

These problems are open for discussions but the impor- 
tant contribution of non thermal UV sources seems to be 
required. Further progress can be achieved with a richer and 
more refined sample of observed metal systems and earlier 
galaxies. 
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